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Abstract 

Exact and approximate expressions are established for dissipation, the power of the 
shear stress at the wall and the boundary layer thickness corresponding to the motion of 
an Oldroyd-B fluid induced by a constantly accelerating plate. The similar expressions for 
Maxwell, Newtonian and second grade fluids, performing the same motion, are obtained 
as limiting cases of our general results. The specific features of the four models are em- 
phasized by means of the asymptotic approximations. 
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1 Introduction 

In recent years, the interest for flows of Newtonian and non-Newtonian fluids has 
considerably increased, while energetic aspects regarding these motions are scarcely met 
in the literature. The first results of this kind seem to be those of Biihler and Zierep [1] 
concerning the Rayleigh-Stokes problem for Newtonian fluids. These results have been 
recently extended by Zierep and Fetecau [2, 3] and Fetecau et al [4] to second grade. 
Maxwell and Oldroyd-B fluids. In their works, the authors also determine the boundary 
layer thickness and use it to obtain a series solution for the velocity field. Similar results 
for the flow induced by a constantly accelerating plate in second grade and Maxwell fluids 
have been also obtained by Fetecau et al in [5, 6]. For this motion, unlike the previous 
flow, the series solutions that have been obtained are completely determined by means of 
the corresponding initial and boundary conditions. 
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The aim of this paper is to extend the results from [5] and [6] to Oldroyd-B fluids. 
More exactly, we intend to present a complete study of the energetic balance for the 
unsteady flow of such a fluid, driven by the transversal motion of a constantly accelerating 
plate. Of special interest are: the Power due to the shear stress at the wall. Dissipation 
and the Changing of the kinetic energy with time. The first term describes the energy 
input that is necessary to keep the medium running. A decisive question is whether these 
terms are larger or smaller than in Maxwell, Second grade and Newtonian fluids. In order 
to see this, besides their exact values, approximative expressions have been established 
for small dimcnsionlcss relaxation and retardation times. Such expressions have been also 
achieved for velocity, shear stress and boundary layer thickness. Finally, in contrast with 
the previous results concerning the flow due to a suddenly moved plate [2-4, 7] and in 
accordance with those obtained in [5] and [6], the series solutions for velocity and shear 
stress have been completely determined by governing equations and the corresponding 
initial and boundary conditions. In the special cases, when the relaxation or retardation 
time tends to zero, all results that have been obtained reduce to those for second grade, 
respectively. Maxwell fluids. If both times are going to zero, the results for Newtonian 
fluids are recovered. 



2 Governing equations 

Recently, Rajagopal and Srinivasa [8] have built and developed a systematic thermo- 
dynamic framework within which models for a variety of rate type viscoelastic fluids can be 
obtained. Among them, the Oldroyd-B fluids that take into account elastic and memory 
effects exhibited by most polymeric and biological liquids have been used quite widely in 
many applications and the results of simulation flt experimental data in a wide range [9] . 
These fluids can describe stress-relaxation, creep and the normal stress differences that 
develop during simple shear flows. They store energy like a linearized elastic solid, their 
dissipation however being due to two dissipativc mechanisms that implies that they arise 
from a mixture of two viscous fluids. The Cauchy stress in an incompressible Oldroyd-B 
fluid, is given by [8-10] 

T = -pI + S, S + A(S-LS-SL^) ==/x[A + A^(A-LA- AL^)], (1) 

where —pi denotes the indeterminate spherical stress, S is the extra-stress tensor, L is the 
velocity gradient, A = 'L+'L^ is the flrst Rivlin-Ericksen tensor, /j is the dynamic viscosity, 
A and A^ are the relaxation and retardation times and the superposed dots indicate the 
material time derivatives. This model includes, as a special cases, the Maxwell and linearly 
viscous fluid models. Let us consider an incompressible Oldroyd-B fluid at rest, occupying 
the space above an infinitely extended plate which is situated in the plane y = of a 
Cartesian coordinate system x, y and z. At time t = 0+, the plate is subject to a motion 
of constant acceleration A in the a;-direction. Owing to the shear the fluid above the plate 
is gradually moved, its velocity being of the form 

V = v{y,t) = u{y,t)i, (2) 

where i denotes the unit vector along the x-coordinate direction. For this flow, the con- 
straint of incompressibility is automatically satisfied. Assuming that the extra-stress S 
depends only of y and t and having in mind the initial condition 

S{y,t)=0, y>0, (3) 
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the fluid being at rest up to the moment t = 0, One obtains (see for instance [11]) 

Sxz — Syz ~ ^yy ~ ^zz ~ 

As regards the shear stress T{y,t) = Sxy{y,t) satisfies the partial differential equation 

(1 + XdtMy, t) = m(1 + Xrdt)dyu{y, t). (4) 

In the absence of pressure gradient in the flow direction and neglecting body forces, 
the balance of linear momentum leads to the meaningful equation 



dyT{y,t) = pdtu{y,t); y,t > 0, 



(5) 



Where p is the constant density of the fluid. Eliminating T{y,t) between Eqs. (4) and (5), 
one attain to the governing equation 



Xdfu{y,t) + dtu{y,t) = 1^(1 + Xrdt)d^u{y,t); y,t>0, 



(6) 



where u = p/p is the kinematic viscosity of the fluid, p being its constant density. 

The exact solutions of the partial differential equations (4) and (6) with the initial 
and boundary conditions 



u{y,0) = dtu{y,0)=0, r(y,0) = for y > 0, 
u{0,t) = At for t>0; u{y,t), dyu{y,t) ^ as 2/ — > 00. 



are given by (cf. with [11], Eqs. (16) and (20)) 



and 



2A f 

u{y,t)=At / 

z^Ti" Jo 

, , 2pA 
r{y,t) = -^ 

TT 

ri,2 = 



^ _ r^r?. exp(ri^) - rir^ exp(r2^) ^ 



r-2 - r-i 



sin(yO 



1 - 



Ti exp(r2t) - rs exp(rit) 



r2 - ri 



cos(y^ 
^2 



(7) 
(8) 

(9) 
(10) 



-(1 + g^^) ± v/(l + a^Y - 4.vXe 
2A ' 



1 - ± v^(l + a^Y - 4^^ 
2A 



and a = vXr- 

By letting or A ^ into Eqs. (9) and (10), ), we attain to the similar solutions 
for Maxwell (see [11], Eqs. (25) and (26)) and second grade (see [12], Eq. (2.5) with 
V{t) = At or [13], Eq. (6) for the velocity field only) fluids. If both A^ and A ^ 0, the 
solutions [11] 

sin(2/0 



UN{y,t) = At / 1 

2^7r Jo V 

2pA 



TN{y,t) = -- 



1 - e--«^* 



e 



(11) 



corresponding to a Newtonian fluid. The last solutions can be also obtained from (9) and 
(10) by making A = A^. This is not a surprise, it is a simple consequence of the Joseph's 
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remark ([14], §2.2). Of course, UM{y,t) and TN{y,t) given by Eqs. (9) can also be written 
in the equivalent classical forms (cf. [65], Eqs. (37) and (38)) 

UN{y,t) = 4:Ati^EYic( a.nd TN{y,t) = -2pAViA iErfcl —^], (12) 

where i"Erfc(a;) = i"-^Erfc(C) d£, are the integrals of the complementary error function 
Erfc(-) and i°Erfc(a;) ^ Erfc(3;). 

The energetic balance for a given volume V, as it results from [1-6], is given by 



di 



Ekin + L + ^ = 0, 



(13) 



where Ekin is the kinetic energy, $ is the dissipation and L is the power of the shear stress 
at the wall. The changing of the kinetic energy with time is given by 



(14) 



where A is the boundary of the flow domain and n is the unit vector normal to A. For 
an infinite volume of rectangular cross-section with x G [0, 1] and z G [0, 1], Eqs. (2), (5), 
(14) imply 

j^Ekin = p J^uiy,tf-^dV = pi f uiy,tf-^dy = lf u^'-I^dy = 



-lu%t)T{Q,t)-l r T{y,t)^^^^dy. 
Jo 



dy 



(15) 



Comparing Eqs. (13) and (15), it results that 



L = L(t) = /u(0,t)r(0,t) = ;w„(i)T„(t) and $ = $(t) = / / T{y,t) ^^^^'*^ dy. (16) 



(17) 



/o dy 

The boundary layer thickness, as it results from [7], is given by 

1 



S = S{t) 



i{y,t)dy 



and represents the thickness of the fluid layer moved with the plate by friction. One 
measure of the boundary layer thickness is the distance from the wall where the velocity 
differs by 1 percent from the external velocity. 



3 Exact expressions for L, $ and 5 

Introducing the general solutions u{y,t) and T{y,t), given by Eqs. (9) and (10), into (16) 
and (17) we find for L, $ and S the exact expressions 



L = 

AplA^ 



2piAH r 

Jo 



Ti exp(r2i) - rs exp(rit) 



1 

Ti exp(r2t) - ra exp(ri<) 
r2 - ri 



(18) 
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r2r3exp(rii) - rir4exp(r2i) 



r-2 - n 



cos(yO 



d^Uy 



and 
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^ _ r2r3 exp(rit) - nr^ exp(r2t) ^' 
r2 - ri 



(19) 



^^U'i.- (20) 



By letting and A ^ or A ^ into Eqs. (18)-(20), we attain to the similar 
expressions (cf. with Eqs. (3.4)-(3.6) from [6]) 



'^2 ' 



JV 



and 



5n = 



Jo 

Jo VJo s . 

sin(y^) 



2 



(1 - e-"'^ 



dy 



dy 



(21) 
(22) 

(23) 



corresponding to a Newtonian fluid performing the same motion. The similar expressions 
for Maxwell or second grade fluids are also obtained as limiting cases of (18)-(20) for A^, 
respectively, A ^ 0. The integrals in Eqs. (21)-(23) can be found and Ln, $jv and 6n 
can be written in simple forms 



-2plA'U -, 



N 



and 



KV2-1) 



plAH\ —, 

IT 



(24) 



in which the boundary layer thickness is independent of the constant acceleration A. 
Indeed, using Eq. (A3) from appendix with y = 0, Eq. (24) i is immediately obtained. In 
order to get the simplified form of <I>Ar, we must use Eqs. (A3), (Al)j^, (A4) and (A5). As 
regards the boundary layer thickness, from (23), (A6) and (A7) one obtains (24)3. 



4 Asymptotic approximations 

4.1 \/t and \/t <C 1 

Let us now consider the case when the dimensionless relaxation and retardation times 
X/t and \r/t = al{vt) are much less then one. In this case the terms containing exp(r2t) 
from Eqs. (9), (10) and (18)-(20) can he neglected, they going to zero faster then (A/i)^. 
. Furthermore, for convenience (in order to get the results for second grade fluids as a 
limiting case) we shall use a/{yt) instead of Xr/t. In these conditions, the next approxi- 
mations 

^/{iT^^^^f^^^i^ = 1 + (a - 2i/A)^2 - 2vX{a - i/A)^"* + • • •, 
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^rit _ „-vtf 



e' - = e (1 + avtC - v^tX^ + •••); 



^•3 



r2 - ri 

t2 1^ I ,,\t2/i ,,v.c2\ 



r-i - ri 



= -e-^*^ [1 + oc{vti' - 1) + v\i\\ - vte) + •••], 



Ar2r3e''i* 



r2 - ri 

are valid for each ^ and t greater than zero. 

Introducing (24)5^6 into (5) and (6) and having Eqs. (9) in mind, we find that 

t) = UN{y, t) A / ^e-'^*^ ^HvOd^ + ■■■ 



and 



T{y,t) = TN{y,t) + ^A /"(I - cosivOd^ + ' 

I" Jo 



(25) 



(26) 



(27) 



Using now Eqs. (A. 2)2 , (A. 8) and (A. 9) from appendix, we can write the simpler forms 



u{y,t) = UN{y,t) - —\ — exp ~ -rz ] ■ T + ^ 



2 y i^n 



and 



where 



r{y,t) = TN{y,t) + 



iiA / t 



2 y VTT 



1 + 



21^1 



Avt j t 

y'\ A 



exp - 



Ai^tJ t 



+ 



(28) 



(29) 



P = Max{a/{vt),\/t}. 



5 Conclusion 

In this paper it is presented a stndy of the energetic balance corresponding to the unsteady 
flow of an Oldroyd-B fluid due to a constantly accelerating plate. Exact and approximative 
expressions are established for dissipation, the power due to the shear stress at the wall 
and the boundary layer thickness. As a conHc;quence, the changing of the kinetic energy 
with time is obtained from the energy balance. In the special cases, when the relaxation 
or retardation time tends to zero, our results are going to the corresponding results for 
second grade and Maxwell fluids. If both times are going to zero, the similar results for 
Newtonian fluids are recovered. 

Finally, by analogy with the Teipel's series expansion [7], series solutions have been 
established both for the velocity filed and the adequate shear stress . Our solutions, un- 
like the series solution of Teipel, are completely determined by means of the governing 
equations and the appropriate boundary conditions. Furthermore, as it was to be ex- 
pected, these solutions are identical to those resulting from exact solutions by means of 
the asymptotic approximations. 

In comparison or by analogy with the unsteady flow induced by a suddenly moved 
plate in Newtonian and non-Newtonian fluids, the next remarks can be made: 
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1. The boundary layer thickness corresponding to these motions is independent of the 
constant velocity U or the constant acceleration A of the plate for all models. 

2. If A > Ar, in comparison with Newtonian fluids, L and (j) increase for the flow in- 
duced by a suddenly moved plate and decrease for the flow produced by a constantly 
accelerating plate, while S decreases for both motions. 

3. The series solutions corresponding to the flow induced by a suddenly moved plate 
contain two free constants, although all conditions have been fulfilled, while the 
series solutions for the flow due to a constantly accelerating plate are completely 
determined by moans of the initial and boundary conditions. 

4. The approximative solutions obtained by means of the series expansions are iden- 
tical to those resulting from the exact solutions using asymptotic approximations. 
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